Measurement of single ir^ production in neutral current neutrino interactions with 
water by a 1.3 GeV wide band muon neutrino beam 
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Neutral current single tt" production induced by neutrinos with a mean energy of 1.3 GeV is 
measured at a 1000 ton water Cherenkov detector as a near detector of the K2K long baseline 
neutrino experiment. The cross section for this process relative to the total charged current cross 
section is measured to be 0.064±0.001 {stat.)±0.007 (sys.). The momentum distribution of produced 
TT^s is measured and is found to be in good agreement with an expectation from the present knowledge 
of the neutrino cross sections. 

PACS numbers: 13.15.+g, 14.60.Lm, 25.30.Pt 



After the discovery of atmospheric neutrino oscillations 
by Super-Kamiokande in 1998 the primal aim of cur- 
rent and future long baseline (LBL) experiments using 
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an accelerator-based neutrino beam is more accurate de- 
termination of oscillation parameters. The uncertainties 
on the knowledge of the neutrino-nucleus cross sections 
and subsequent nuclear effects in the GeV neutrino en- 
ergy region could become a severe limitation in future 
oscillation studies. For reducing these systematic uncer- 
tainties, near detectors of LBL experiments can provide 
neutrino interaction data with much higher statistics and 



2 



better quality than existing data. 

A single tt" event is a good signature of neutral current 
(NC) neutrino interactions in the GeV region. Especially 
in a water Cherenkov detector, a decay of the tt" can 
be clearly identified as two electromagnetic-showering 
Cherenkov rings. The single tt'^ production rate by atmo- 
spheric neutrinos could be usable to distinguish between 
the h'f^ <-^- i^r and <-> Vs oscillation hypotheses. The NC 
rate is attenuated in the case of transitions of i/^ 's into 
sterile neutrinos, while it does not change in the <-!■ 
scenario. An understanding of single tt*' production is 
also very important for a search for electron neutrino ap- 
pearance in LBL experiments with a water Cherenkov 
detector, because the most serious background to single- 
ring signals is a single 7r° event with only one ring re- 
constructed due to highly asymmetric energies or small 
opening angle of two 7-rays in the decay • 

These single tt'^s in the GeV neutrino energy region 
are mainly produced via the A resonance as v + N 
1/ + A, A ^ iV' + tt", where N and N' are nucleons. 
Because of nuclear effects such as Fermi motion, Pauli 
blocking, nuclear potential and final state interactions, A 
production and its decay could be different from a sim- 
ple picture of neutrino- nucleon interactions. In addition, 
final state interactions of nucleons and mesons during 
a traversal of nuclear matter could largely modify the 
number, momenta, directions and charge states of pro- 
duced particles. Though there exist several theoretical 
approaches for modeling these processes, their uncertain- 
ties are still large. There exist very little experimental 
data for NC single 7r° production and no reports on mea- 
surements with a water target, which is the target matter 
of a far detector in some of the LBL experiments 0, 01 ■ 
It is clear that a good knowledge of NC single 7r° pro- 
duction cross section and 7r° momentum distribution is 
required for the above studies. In this letter, we present 
the first high statistics measurement of " NCItt" " in- 
teractions in water defined as a neutral current neutrino 
interaction where just a single tt^ and no other mesons 
are emitted in the final state from a nucleus. 

The KEK to Kamioka long baseline neutrino oscilla- 
tion experiment (K2K) uses an almost pure muon neu- 
trino beam (98.2% f^, 1.3% i/e and 0.5% j7^) with a 
mean neutrino energy of 1.3 GeV. This intense wide-band 
beam is produced by the KEK proton synchrotron (KEK- 
PS). Protons with 12 GeV kinetic energy extracted from 
the KEK-PS are bent toward the far detector, Super- 
Kamiokande, located 250 km away from KEK and inter- 
act with an aluminum target. Positively charged sec- 
ondary particles, mainly tt+'s, are focused with a pair 
of magnetic horns and then decay to produce a neutrino 
beam. Figure ^ shows the energy spectrum of the K2K 
neutrino beam at 300 m downstream from the target in 
the near site with a 10^° protons on target (POT) expo- 
sure predicted by a beam simulation 01 . The beam sim- 
ulation is validated by a pion monitor (PIMON), which 
measures the momentum and divergence of pions just be- 
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FIG. 1: The energy spectrum of the K2K neutrino beam at 
300 m downstream from the target in the near site with a 
10^" protons on target exposure predicted by a neutrino beam 
simulation. The spectrum is averaged within 2 m from the 
beam center. 



hind the second horn (3| . The flux shape is also measured 
through neutrino interactions by K2K near detectors |^ . 
The absolute flux estimation, however, is not easy due 
to uncertainties of the primary proton beam absolute in- 
tensity, the proton beam profile and the proton targeting 
efficiency. Instead of deriving the absolute cross section 
for NCItt^ interactions, we measure the relative NCItt" 
cross section to the total charged current (CC) cross sec- 
tion, since the CC is a good experimental signature. 

As one of the near detectors for K2K, a 1,000 ton water 
Cherenkov detector (Ikt) is located about 300 m down- 
stream of the pion production target. The Ikt detector 
is a miniature of Super-Kamiokande using the same in- 
teraction target matter and instrumentation. The inner 
volume of the Ikt detector is a cylinder with 8.6 m diam- 
eter and a height of 8.6 m. This volume is viewed by 680 
photomultiplier tubes (PMTs) of 50 cm diameter, facing 
inward to detect Cherenkov light from neutrino events. 
The PMTs and their arrangement are identical to those 
of Super-Kamiokande, giving a 40 % photocathode cover- 
age. The primal role of the Ikt detector is to measure the 
i^^ interaction rate and the energy spectrum. The Ikt 
detector also provides a good measurement of high statis- 
tics neutrino-water interactions. The analysis in this let- 
ter is based on the Ikt data taken between January 2000 
and July 2001, corresponding to 3.2 x lO^^ POT. 

The data acquisition (DAQ) system of the Ikt detec- 
tor is also similar to that of Super-Kamiokande. A signal 
from each PMT is processed by custom electronics mod- 
ules called ATMs, which are developed for the Super- 
Kamiokande experiment and are used for recording dig- 
itized charge and timing informations of each PMT hit 
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over a threshold of about 1/4 photoelectrons '6!|. The 
DAQ trigger threshold is about 40 hits of PMTs within 
a 200nsec window in a 1.2/xsec beam spill gate, which is 
roughly equivalent to a signal of a 6MeV electron. The 
analog sum of aU 680 PMTs' signals (PMTSUM) is also 
recorded for every beam spill by a 500 MHz FADC to 
identify multi interactions in a spill gate. We determine 
the number of interactions in each spill by counting the 
peaks in PMTSUM greater than a threshold equivalent 
to a 100 MeV electron signal. 

Physical parameters of an event in the Ikt detector 
such as the vertex position, the number of Cherenkov 
rings, particle types and momenta are determined by 
the same algorithms as used in Super-Kamiokande [3, Q ■ 
First, the vertex position of an event is determined by 
timing information of PMT hits. With knowledge of the 
vertex position, the number of Cherenkov rings and their 
directions are determined by a maximum-likelihood pro- 
cedure. Each ring is then classified as e-like representing 
a showering particle (e^,7) or /x-like representing a non- 
showering particle (/i^,7r^) using its ring pattern. On 
the basis of this particle type information, the vertex 
position of a single-ring event is further refined. The mo- 
mentum corresponding to each ring is determined from 
the intensity of Cherenkov light. 

In this analysis, neutrino interaction candidates are 
selected by the following requirements : (i) An event is 
triggered within a 1.2/isec beam spill gate, (ii) There 
is no detector activity within a 1.2 ^sec window preced- 
ing a beam spill, (iii) Only a single event is observed by 
the PMTSUM peak search in that spiU. (iv) The recon- 
structed vertex position is in the 25 ton fiducial volume 
defined as a 4m diameter, 2m long cylinder along the 
beam axis, (v) The visible energy is larger than 30 MeV. 
A total of 60,545 events remain after cuts (i) to (v). 

Single tt'^ events are extracted from the sample of fully 
contained (FC) neutrino events, which deposit all of their 
Cherenkov light inside the inner detector. An exiting par- 
ticle deposits a large amount of energy at the PMT of the 
exiting position. Therefore, the FC events are selected by 
requiring the maximum number of photoelectrons on a 
single PMT at the exit direction of the most energetic 
particle to be less than 200. The events with the maxi- 
mum number of photoelectrons greater than 200 are iden- 
tified as a partially contained (PC) event. Because tt^s 
mostly decay into two 7-rays, the following criteria are 
further applied to the FC sample in order to select sin- 
gle TT*^ events ; (1) The number of reconstructed rings is 
two. (2) Both rings have a showering type (e-like) parti- 
cle identification. (3) The invariant mass is in the range 
of 85 — 215MeV/c^. Figure 13 shows the invariant mass 
distribution of the events which satisfy cuts (1) and (2). 
A 7r° mass peak is clearly visible. The peaks for both the 
observed data and the neutrino Monte Carlo events are 
slightly shifted towards higher values compared to the 
nominal value of the 7r° mass, 135MeV/c^. This shift 
is due to a vertex reconstruction bias of a few dozen cm 
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FIG. 2: The invariant mass of two e-like ring events for the 
experimental data (black dots) and the neutrino Monte Carlo 
simulation (box histogram). The Monte Carlo histogram is 
normalized to the area of the data histogram. The error bars 
are statistical only. The black portion and the hatched por- 
tion in the Monte Carlo histogram show the non-Tr'' com- 
ponent and non-NClTr" component, respectively. A pair of 
arrows shows the invariant mass cut (3) (see text). 





Data 


NCItt" efficiency 


FC 


45317 


97% 


Two rings 


11117 


57% 


Both e-hke 


3150 


48% 


Invariant mass 


2496 


47% 



TABLE I: The number of events after each selection to make 
the single n° sample in Ikt data. The Monte Carlo efficiencies 
are calculated for NCItt'^ interactions whose real vertex is in 
the fiducial volume. 



and an energy deposit by de-excitation 7-rays from an 
oxygen nucleus. The difference in peak position between 
the data distribution and the Monte Carlo prediction is 
within our quoted I3 % systematic uncertainty for the 
absolute energy scale. A total of 2,496 single 7r° events 
are collected by these criteria as shown in Table 

Figure |Sl-(a) shows the reconstructed tt" momentum 
distribution for the single tt'^ sample. The momentum 
resolution for 200MeV/c 7r"s in NCItt" events is es- 
timated to be about 15MeV/c. The single tt" sam- 
ple contains a background of non-NClTr" events shown 
as a hatched histogram in Figures |21 and |31-(a). The 
background contamination is estimated by the neutrino 
Monte Carlo simulation and is subtracted. 

The neutrino interactions with water are simulated by 
the NEUT program library. A detailed description of 
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FIG. 3: (a) The reconstructed tt" momentum distribution for 
the single ■k'^ sample, comparing Ikt data (black dots) and 
the neutrino Monte Carlo simulation (box histogram). The 
Monte Carlo histogram is normalized to the area of the data 
histogram. The error bars are statistical only. The black por- 
tion and the hatched portion in the Monte Carlo histogram 
show the non-7r° component and non-NClTr" component, re- 
spectively, (b) The detection efficiency for NCItt" interactions 
as a function of real 7r° momentum. The highest momentum 
bin in each figure integrates the events above 800MeV/c. 



NEUT can be found in Ref . . Quasi-elastic scatte ring 
is simulated based on the Llewellyn Smith's model [lO|. 
The Rein and Sehgal's model [Ij is used to simulate 
single meson production mediated by a baryon reso- 
nance. The decay kinematics of the A(1232) resonance 
is also determined by the Rein and Sehgal's method. 
For the decays of the other resonances, the generated 
mesons are assumed to be emitted isotropically in the 
resonance rest frame. A twenty percent suppression of 
pion production is adopted for simulating pion-less A 
decay, where the event contains only a lepton and a nu- 
cleon in the final state The axial vector mass. Ma, 
is set to l.lGeV/c^ for both the quasi-elastic scatter- 
ing and the single meson production models. Coherent 
pion production is simulated using the Rein and Seh- 
gal's model JJ3J with modified cross section according to 
a description by Marteau et al. For deep inelastic 

scattering, GRV94 nucleon structure function with 
a correction by Bodek and Yang is used to calcu- 
late the cross section. In order to generate the final state 
hadrons, PYTHIA/JETSET ^ package is used for the 
hadronic invariant mass, W, larger than 2.0GeV/c^. A 
custom-made program [18| is used for W in the range 
of 1.3 - 2.0GeV/c2. Nuclear effects for v-'^^O scatter- 



ing are also taken into account. Fermi motion of nucle- 
ons is simulated using the relativistic Fermi gas model. 
Pauli blocking effect is considered in the simulation of 
quasi-elastic scattering and single meson production. Pi- 
ous generated in ^^O are tracked taking into account 
their inelastic scattering, charge exchange and absorp- 
tion. The pion generation point in the nucleus is set 
according to Woods and Saxon's nucleon density distri- 
bution jl^l . The mean free path of each pion interaction 
is calculated using the model by Salcedo et al. The 
direction and momentum of pious after inelastic scatter- 
ing or charge exchange are determined based on the re- 
sults of a phase shift analysis obtained from w — N scat- 
tering experiments Emission of de-excitation low 
energy 7-rays from a hole state of residual nuclei is also 
taken into account ^22]. Outside the nucleus, particles 
are tracked using GEANT and CALOR packages except 
pions with momenta smaller than 500 MeVV c, which are 
tracked with a custom-made program jlSj in order to 
consider the A resonance correctly. 

The NCl7r° fraction in the single 7r° sample, predicted 
by the neutrino Monte Carlo simulation, is estimated to 
be 71 % (52 % from single tt" production via a resonance 
decay, 3 % from single tt^ production via a resonance de- 
cay and subsequent charge exchange into tt^ in the tar- 
get nucleus, 10% from coherent tt" production and 4% 
from neutrino deep inelastic scattering where the rest 
of mesons are absorbed inside the nucleus). The non- 
NCItt" background, which accounts for 29% of the sin- 
gle tt" sample as shown in Figures El and |31 contains NC 
tt'^ production where outgoing mesons except a single tt*' 
have low momenta (7%), CC tt^ production where ac- 
companying muon and mesons have low momenta (9%), 
TT^ production by a recoil nucleon or mesons outside the 
target nucleus (10%), and non-7r" background due to 
mis-reconstruction (3%). 

The NCItt^ fraction is estimated for each of the ten 
tt'^ momentum bins shown in Figure 13- (a). The number 
of observed single 7r° events in each 7r° momentum bin 
is multiplied by the NCItt" fraction for the correspond- 
ing momentum bin. The systematic errors attributed 
to this background subtraction procedure are estimated 
by assuming uncertainties in neutrino interaction models 
as listed in Table IIll-(A). The Ma value in the models 
for quasi-elastic scattering and single meson production 
is varied by ±10% from a central value of l.lGeV/c'^. 
The error due to the uncertainty on coherent pion pro- 
duction cross section is estimated by removing the cross 
section modification by Marteau et al.. For deep inelastic 
scattering, the effect of the correction on nucleon struc- 
ture function according to Bodek and Yang's description 
is evaluated. The systematic errors due to nuclear effect 
uncertainties are evaluated by varying the probabilities of 
pion absorption and inelastic scattering (including charge 
exchange) in ^^O nucleus independently by ±30%. 

After estimating the number of NCItt'' interactions 
in each tt'^ momentum bin, we apply a fiducial volume 
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Sources Errors (%) 

(A) Systematic uncertainties in background subtraction 

Ma in quasi-elastic and single meson (±10%) 0.2 

Quasi-elastic scattering (total cross section, ±10 %) 0.0 

Single meson production (total cross section, ±10%) 0.9 

Coherent pion production (model dependence) 1.6 

Deep inelastic scattering (model dependence) 5.1 

Deep inelastic scattering (total cross section, ±5 %) 0.5 

NC/CC ratio (±20 %) 3.2 

Nuclear effects for pions in ^^O (absorption, ±30%) 1.5 

Nuclear effects for pions in ^"O (inelastic scattering, ±30%) 0.7 
Pion production outside the target nucleus (total cross section, ±20 %) 2.3 

(B) Systematic uncertainties in fiducial volume correction 

Fiducial cut 1.6 

(C) Systematic uncertainties in efficiency correction 

Ring counting 5.4 

Particle identification 4.2 

Energy scale 0.3 



TABLE 11: Summary of the systematic errors on the measurement of the number of NCItt interactions. 



correction. The number of NCItt^ interactions in the 
Monte Carlo single tt'^ sample increases by about 2 % if 
real vertices are used instead of reconstructed vertices in 
the fiducial volume selection. The data are then multi- 
plied by 1.02 in order to derive the number of NCItt" 
interactions in the true 25 ton fiducial volume before the 
efficiency correction is applied. The systematic error in 
the fiducial volume correction is about 2 % as shown in 
Table m-(B), which is estimated from the difference in 
the reconstructed vertex distribution of single tt" events 
for Ikt data and the neutrino Monte Carlo simulation. 

Finally, the number of NCItt" interactions in the true 
25 ton fiducial volume is corrected for the detection effi- 
ciency in a bin by bin manner. Figure 13- (b) shows the 
detection efficiency for NCfTr" interactions as a function 
of real tt" momentum. The inefficiency for higher mo- 
mentum Tr"s is primarily due to reconstruction of only 
one ring for the tt" decay with highly asymmetric ener- 
gies or small opening angle of two 7-rays. The overall 
detection efficiency for NCf tt" interactions, estimated by 
the neutrino Monte Carlo simulation, is 47% as shown in 
Tabled The systematic errors from the efficiency correc- 
tion are due to uncertainties of reconstruction algorithms 
such as ring counting, particle identification and energy 
scale as listed in Table llll-(C'). 

By a series of corrections for the background subtrac- 
tion, the true/reconstruction fiducial difference and the 
detection efficiency, the true number of NCItt" interac- 
tions in our data is measured to be (3.61 ± 0.07 ± 0.36) x 
lO'^ in the 25 ton fiducial volume. Figure 0] shows the 
measured momentum distribution of NCItt" after all cor- 
rections, compared with the distribution predicted by the 
neutrino Monte Carlo simulation. The Monte Carlo his- 
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FIG. 4: The momentum distribution of NCItt" events in the 
25 ton fiducial volume (black dots) . The inner and outer error 
bars attached to data points show statistical errors and total 
errors including systematic errors, respectively. The distri- 
bution predicted by the neutrino Monte Carlo simulation is 
also shown as a box histogram for comparison. The size of in- 
ner boxes represents Monte Carlo statistical errors. The size 
of outer boxes represents the uncertainty of the distribution 
shape due to neutrino interaction model ambiguity. 



togram is normalized by the number of total neutrino 
events in the fiducial volume, which are selected by cuts 
(i) to (v) as previously described. The size of outer 
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boxes for the Monte Carlo histogram represents the un- 
certainty of the distribution shape of NCItt*^ momentum 
due to neutrino interaction model ambiguity, where the 
largest source is nuclear effects for pions in -'^^O. The mea- 
sured distribution is in reasonably good agreement on the 
Monte Carlo estimation. 

As previously described, we use CC interactions for 
normalization in order to derive the relative cross section 
for NCl7r° interactions. To make a CC enriched sample, 
FC /x-like events and PC events are selected from the Ikt 
data during the same period. By using the CC enriched 
sample as a normalization, the uncertainty of the neu- 
trino energy spectrum is almost canceled in the mea- 
surement since the expected mean energy of neutrinos 
producing the CC sample, 1.45 GeV, is almost same as 
that of neutrinos producing the 7r° sample, 1.50 GeV. The 
sample consists of 22,612 FC single-ring /i-like events, 
12,386 FC multi-ring events with the most energetic ring 
identified as /x-like and 15,228 PC events, resulting in a 
total of 50,226 events. The t'^CC fraction in this sam- 
ple is estimated to be 96 % by the neutrino Monte Carlo 
simulation (96.5% for the FC single-ring ^-like sample, 
91.2% for the FC multi-ring /u-like sample and 98.5% 
for the PC sample). The rest 4 % of the sample is mostly 
composed of NC interactions accompanying an outgo- 
ing charged pion above its Cherenkov threshold. The 
fiducial volume correction factor is estimated to be 1.02. 
The detection efficiency for i/^CC interactions by this se- 
lection is estimated to be 85%. The inefficiency mainly 
comes from mis-identification of the ring type in multi- 
ring events and ^ 100 MeV visible energy threshold by 
peak counting of the PMTSUM signal. 

By applying non-i'^CC backgound subtraction, fidu- 
cial volume correction and detection efficiency correction, 
the number of z^^CC neutrino interactions during the 
analysis period is measured to be (5.65±0.03±0.26) x 10^ 
in the 25 ton fiducial volume. The estimated systematic 



errors are 4 % from the uncertainty of vertex reconstruc- 
tion, 1 % from the uncertainty of neutrino interaction 
models, 1 % from the uncertainty of particle identification 
and 1 % from the uncertainty of absolute energy scale. 

By taking the ratio, the relative cross section for 
NCItt^ interactions to the total f^CC cross section 
is measured to be 0.064 ± 0.001 {stat.) ± 0.007 (sys.). 
Our neutrino interaction models predict the ratio to 
be 0.065, which results in good agreement. For refer- 
ence, the total J^^CC cross section is calculated to be 
1.1 X 10"'^* cm^/nucleon in the neutrino Monte Carlo sim- 
ulation by averaging over the K2K neutrino beam energy. 

In summary, we have measured the rate and the tt^ 
momentum distribution of NC single 7r° production by 
neutrinos with a mean energy of 1.3 GeV at the K2K 
Ikt water Cherenkov detector. The cross section ra- 
tio to total i^^CC neutrino interaction is obtained as 
0.064 ± 0.001 {stat.) ± 0.007 (sys.), showing good agree- 
ment with the prediction by our neutrino interaction 
models. This measurement provides essential informa- 
tion for an understanding of 7r° production in water at 
the 1 GeV neutrino energy region, which is relevant for 
present and future neutrino oscillation experiments. 
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